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Abstract 

The aim of this paper is to provide new results and insights for a screw dislocation 
in functionally graded media within the gauge theory of dislocations. We present 
the equations of motion for dislocations in inhomogeneous media. We specify the 
equations of motion for a screw dislocation in a functionally graded material. The 
material properties are assumed to vary exponentially along the x and y-directions. 
In the present work we give the analytical gauge field theoretic solution to the prob- 
lem of a screw dislocation in inhomogeneous media. Using the dislocation gauge ap- 
proach, rigorous analytical expressions for the elastic distortions, the force stresses, 
the dislocation density and the pseudomoment stresses are obtained depending on 
the moduli of gradation and an effective intrinsic length scale characteristic for the 
functionally graded material under consideration. 
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1 Introduction 



Nonhomogeneous media, multilayered structures and functionally graded materials (FGMs) 
are of consid erable technical and engineering importance as well as of mathematical inter- 
est (see, e.g.. lErdoganI ( 1l995l )). Generally, such materials refer to heterogenous composite 



materials in which the material moduli vary smoothly and continously from point to point. 
Typical examples of FGMs are ceramic / cer amic and metal/ceramic s yste ms. Comp ressive 
reviews on aspects of FGMs can be found in lMarkworth et al.l (119951 ) and lErdogaru (11995 ). 
However, only few investigations have been made to assess the role and the importance 
of dislocations in FGMs. For the first time, iBarnettI (119721 ) found the displacement and 
stress fields of a screw dislocation in an isotropic medium which is arbitrary graded in 
the a;-direction. One important example for a FGM are exponentially graded materials. 
For such media the elasticity moduli are exponential functions depending on t he space 



coordi nates and the new material parameters describing the gradation (see, e.g., lErdogan 
( 119951 )). The Gree n function for a tw o dimensional exponentially graded elastic medium 
ha s been found by I Chan et al.l (120041 ) and the three-dimensional Green function is given 
bv iMartin et al.l (120021 ). 

Classical continuum theories possess no intrinsic length scales and, therefore, they are 
scale-free continuum theories which are not able to describe size effects being of impor- 
tance at micro- and nanoscales and near defects. In order to describe such effects general- 
ized continuum theories are needed. Such generalized continuura theories are, for instance, 
strain gradient theories (Mindlin, 1964 ;^ Altan and AifantisI, 199?!; Lazar and Maugin . 2005 ) 
and the dislocation gauge theory (Kadic and Edelen . 19831 : Edelen and Lagoudad . 19881 : 



Lazarl . 120021 : iLazar and Anastassiadisl . 120091 ) which enrich the classical theories with ad- 
ditional materi al lengths. Using str ain gr adient elasti c ity, c racks in FGMs have been 



investigated by iPaulino et al.l ( 120031 ) and IChan et al.l (120081 ). For the first time ever 



the solut ion of a screw dislocation in FGMs, using the theory of gradient elasticity, was 
given by iLazarl ( 120071 ). Gradient elasticity is a theory which is very popular in engineer- 
ing sciences. But not everything is well-understood about the physical importance of 
the higher-order stresses (hyperstresses) . Gradient elasticity should be understood as an 
effective theory since it was invented for compatible deformations. On the other hand, a 
more physically motivated 'gradient-like' theory is the so-called dislocation gauge theory. 
In the dislocation gauge theory both the incompatible elastic distortion tensor and the 
dislocation density tensor are physical state quantities giving contribution to the elastic 
stored energy density. Also for nonhomogeneous media these physical state quantities 
are gauge-invariant. In this approach, the higher-order stress is realized as pseudomo- 
ment stress which is the response to the presence of dislocations. This pseudomoment 
stress is related to the mome r it str es s known from g eneralized elasticities like Cosserat 
elasticity (see, e.g., iNowackil (119861 ): lEringenI ( 1l999l )). Thus, the force stress tensor is 



asymmetric like in Cosserat elasticity. In gen eral, in gradien t elasticity the hyperstress 
(double stress) pos sesses five material modu li ( iMindlinl . Il964j ). In a simplified version of 
gradient elasticity ( 1 Alt an and AifantisI . Il997l ) under ad-hoc assumptions the double stress 
can be simplified. In the dislocation gauge theory such assumptions are not necessary. 
For that reason we use the gauge theory of dislocations in this paper. 

In this paper we study a screw dislocation in FGM in the framework of the translation 
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gauge theory of dislocations. In Section 2 we start with the general framework of the 
translational gauge theory of dislocations and we derive the isotropic three-dimensional 
field equations of dislocations for nonhomogeneous media as well as for exponentially 
graded materials. In Section 3 we specialize in the anti-plane problem and we give the 
gauge-theoretical solution of a screw dislocation in exponentially graded materials. We 
find rigorous and 'simple' analytical solutions for the force stresses, elastic distortion, 
torsion and pseudomoment stresses produced by a screw dislocation in FGMs. Our hope 
is that the gauge solution will provide new physical insight to possible improvements of 
designing FGMs. 



2 Gauge theory of dislocations 

This section introduces the notation and constitutive e quations of the translational gauge 



theory of dislocations (ILazar and Anastassiadid . 120091 ). which will be used to investigate 



a screw dislocation in FGMs. In the dislocation gauge theory, the elastic distortion tensor 
f3ij is incompatible due to the occurrence of a so-called translational gauge field (pij 

Pij = + 4>ij , (1) 

where Ui is the displacement vector. It is important to note that Ui and (pij are just 
the canonical field quantities and they are not unique and not physical state quantities. 
In presence of dislocations, Ui and (pij are discontinuous (or multivalued) fields. The 
gauge field (pij may be identified with the negative plastic distortion. The appearance 
of incompatibility in form of the translation gauge field gives rise to an incompatibility 
tensor in the framework of the translation gauge theory of dislocations which is called the 
dislocation density tensor or torsion tensor. Accordingly, the dislocation density tensor is 
given in terms of the gauge field 

Tijk 'PikJ (pij,k ; -^ij'fc -^ifcj (2) 

and, alternatively, in terms of the elastic distortion 

Tijk f^ik,j t^ij,k ■ (3) 

In the present framework, the elastic distortion tensor and the dislocation density 
tensor T^fc are the physical state quantities. Due to the form of the dislocation density 
tensor ([2]) and ([3]) it fulfills the translational Bianchi identity 

^jkiTijk,i = . (4) 

Eq. (jlj) is the well-known conservation law of dislocations and states that a dislocation 
line cannot end inside the body. 

In the static case of the (linear) translation gauge theory of dislocations, the La- 
grangian density is of the bilinear form 

C= -W = — aij(3ij - ^ HijkTijk , (5) 
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where W denotes the stored energy density. The canonical conjugate quantities (response 
quantities) are defined by 

dC dC 

(^ij • ~aa ' Hijk ■ 2 , Hijj^ Hi]^j , (6) 



ijk 



where a^j and Hijj^ are the force stress tensor and the pseudomoment stress tensor, re- 
spectively. Here, the force stress tensor a^j is asymmetric. The moment stress ten- 
sor Tjjk = —Tjik can be obta i ned f r om the pseudora o ment stress tensor: Tijk = —H^ij-^k 



see 



Lazar and Anastassiadis ( 2009 ): Lazar and Hehl ( 2010l )). They have the dimensions 



SI SI 

[aij] = force/ (length)^ = Pa and [Hijk] = force/length = N/m. 

The Euler-Lagrange equations with respect to the canonical field variables derived 
from the total Lagrangian density C = C{f3ij, Tijk) are given by 



(7) 



We may add to £ a so-called null Lagrangian, Cn = (^ijl^ij, with the 'background' stress 
afj satisfying: cr°- = 0. In terms of the canonical conjugate quantities ([6]), Eqs. (JTj) and 
(E!) take the form 

'^ijj — , (force balance of elasticity) , (9) 
Hijk,k + o'ij = o'ij , (stress balance of dislocations) . (10) 

The force stresses are the sources of the pseudomoment stress. In fact, one can see in 
Eq. ( fTOl) that the source of the pseudomoment st r ess te nsor is an effective force stress 



tensor {aij — cr° ) (see, e.g., lEdelen and LagoudasI ( 119881 )). This effective stress tensor 



which is the difference between the gauge-theoretical stress field and the background 
stress field, drives the dislocation fields. 

The linear, isotropic constitutive relations for the force stress and the pseudomoment 
stress are 

aij = \5ijl3kk + yu(/3ij + (3ji) + i{Pij - Pji) , (11) 
Hijk = ciTijk + C2(Tjki — Tkji) + cz{5ijTiik — SikTiij) . (12) 

Here /x, A,7 are the elastic stiffness moduli and ci,C2,C3 denote the resistivity moduli 
associated with dislocations. The six material moduli have the dimensions: [/x, A,7] = 

SI SI 

force/ (length)^ = Pa and [ci, C2, C3] = force = N. Here, the material moduli are nonhomo- 
geneous that means they depend on the space coordinates: = fj.{x), X = X{x), 7 = 7(0;), 
ci = ci{x), C2 = C2{x) and C3 = c^i^x). The constitutive relations for nonhomogeneous 
media ( |TT|) and ( |T2|) have the same formal form as the relations for homogeneous media. 
The only difference is that the constitutive moduli are not constant. 

The requirement of non-negativity of the stored energy (material stability) W > 
leads to the conditions of semi-positiveness of the constitutive moduli. Particularly, the 
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constitutive moduli have to fulfill the following conditions (ILazar and Anastassiadis 



20091) 



/i>0, 7>0, 3A + 2^>0, (13) 

ci - C2 > , Ci + 2C2 > , Ci - C2 + 2C3 > . (14) 

If we substitute the constitutive equations (ITT]) and ( [T2|) into the Euler-Lagrange 
equations (Q and (ITO!) and use the definition ([3]), we obtain the explicite form of the 
Euler-Lagrange equations for a nonhomogeneous material 

\(3jj^i + + + (/i - + + (/i + 7),jAj + (/U - 7),j/3ii = , (15) 

Cl{f3ikjk — Pij,kk) + C2{(3ji^kk — Pjk,ik + Pkj,ik — Pki,jk) + C3 [5jj(Afc,iA: " A«,fcfc) + Pkk,ji — Pkj,ki\ 
+ Ci^k{f^ik,j — f^ij,k) + C2,k{f3ji,k — f^jk,i + /3fcj,i — f^kij) + Cs^fc 5ij{l^lk,l — (^ll,k) + C^,i{f^kk,j — (^kj,k) 

+ A 5.,(/3,, - /3°,) + (/i + 7)(A, - /3° ) + (/i - - /3° ) = , (16) 

where denotes the 'classical' elastic distortion tensor. It can be seen that the space- 
dependence of the material moduli influences the field equations (|T5l) and f|T6l) due to 
gradients of the material moduli. Eq. f lTSj) is a partial differential equation of first order 
for the elastic distortion as well as for the constitutive moduli. On the other hand, Eq. flT6|) 
contains differential operators of second order acting on the elastic distortion and 'mixed' 
gradients of the constitutive moduli and the elastic distortions. 

Rather than looking for the general inhomogeneous behaviour of the constitutive mod- 
uli, we specify here on an inhomogeneous material, which can find application to function- 
ally graded materials. We assume that the material properties very in a simple, explicite 
manner. Here, we consider exponential variations of the six constitutive moduli in the 
following manner 

— Q2{a-ix+a2y+aiz) ^ _ ^0 ^2{aix+a2y+mz) ^ _ ^0 ^2{a-ix+a2y+as,z) 
= (P^Q^i"-^^+"-^y+°-'i^) — ^Q^{a.ix+a2y+a'iz) _ ^0 g2(aia;+a2y+a32) ^j^y-j 

where A°, 7°, c?, c^, C3, oi, 02 and 03 are the material parameters of the FGM. The 
moduli ai, 02 and 03 are the characteristic parameters for exponentially gradation in x, 
y and 2;-direction and they have the dimensions: [01,02,03] = 1 /length. The parameters 
Oi, 02 and 03 determine the gradation of all six constitutive moduli A, /x, 7, Ci, C2, C3. For 
such an exponentially FGM, the equations ffT^ and flTB]) reduce to 

A(V, + 2a,)/3,-, + (^ + 7)(V, + 2a,) + (/i - 7)(V, + 2a,)/3ji = , (18) 
ci(Vfc + 2afc)(Vj/3ifc - VfcAj) + C2(Vfc + 2ak){Vkl^ji - ^if^jk + ^il^kj - ^jf^ki) 
+ C3 6ij{Vk + 2ak){Vi(3ik - VkPu) + C3(Vi + 2ai){Vjl3kk - VkPkj) 
+ A 6,0kk - Plk) + (/i + 7)(A, - /3° ) + (/i - 7)(/3,. - /3° ) = , (19) 

where Vj denotes the three-dimensional Nabla operator. Eqs. flTSl) and f[T^ are the three- 
dimensional governing equations of dislocations in exponentially graded materials using 
the dislocation gauge theory. 
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3 Gauge solution of a screw dislocation in function- 
ally graded materials 

We now derive the equations of motion for a screw dislocation in FGMs. We consider 
an infinitely long screw dislocation parallel to the 2;- axis with the Burgers vector b = b^. 
The symmetry of such a straight screw dislocation le aves only the f ollowing non- vanishing 



components of the elastic distortion tensor (see, e.g., IdeWitI ( 119731 )): Pzx, Pzy, and for the 
dislocation density tensor: Tzxy 

For the anti-plane problem of a screw dislocation in nonhomogeneous materials, the 
equilibrium condition f lT^ reduces to 

P.,, = -^-^^-^ , J=x,y. (20) 

Using equation (1201) . we obtain from f ll6l) the following equations for the components Pzx 
and (3zy of a screw dislocation in nonhomogeneous materials 



ClPzx^jj ~l~ Ci j /3zx,j [ Cl 7 Ci 

\ + 7 

C2Pzx,jj + C2.jPzx.j — (c2,j — C2 — — ■ ^— 

y ,j + 7 

r B -^r B -fr-r + 



Pzj,x 


- {lJ' + l){(3zx 


- Pi) 


= 0, 


(21) 


Pzj,x 


+ (/W - l)Wzx 


- PL) 


= 0, 


(22) 


Pzj,y 


- {f^ + l){f3zy 


- O 


= 0, 


(23) 


f^zj,y 


+ {f^-l){Pzy 


- P%) 


= 0. 


(24) 



/i + 7 
/ -|- ^) j 

C2Pzy,jj + C2,jPzy,j " [C2,j - C2 ^ _^ ^ ' 

For the two-dimensional anti-plane problem of FGMs, we assume that the constitutive 
moduli are exponentially graded, 

^ = / Q^ia^x+a^y) ^ ^ _ ^0 ^2ia,x+a2y) ^ Ci = C? e^^'^i^+^^^z) ^ C2 = 4 6^^"'''+''^^'^ , (25) 

where 7°, c?, c^, ai and 02 are the material constants of the anti-plane problem. The 
moduli ai and 02 are the characteristic parameters of exponentially gradation in x and 
y-direction, respectively. If we substitute ( |25l) into (I2T1)-( |M1) . the governing equations are 



ci(A + 2a ■V)Pzx - (/i + l){Pzx - PL) = , (26) 

C2(A + 2a ■V)Pzx + (/i - l){Pzx - PL) = , (27) 

Ci(A + 2a ■V)Pzy - (/i + l){Pzy - P%) = , (28) 

C2(A + 2a ■V)Pzy + (/i - l){Pzy " P%) = , (29) 

where A and V denote the two-dimensional Laplacian and the Nabla operator, respec- 
tively. It can be seen in fl2^ -( l2^ that we have four equations for only two components 
Pzx and Pzy. Thus, it is clear that not all four equations are independent. In 

fact, if we compare f l26|) with fl27|) as well as f l28|) with fl29|) . we find the relation between 
the material moduli 

^ = . (30) 

/i + 7 /i-7 
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In addition, we may define the following intrinsic gauge length of the anti-plane problem 
of a screw dislocation 



Cl 



+ 7 V (/^ + 7)° 



(31) 



Evidently, the length scale I is con stant for FGM s using (l25l). This fact is in agreement 
with gradient elasticity for FGMs ( iPaulino et al.l . l2003l : IChan et al.l . l2006l . I2OO8I : iLazar . 



20071 ) where the characteristic gradient length is also constant. Only the constitutive 



moduli entering the constitutive relations are nonconstant. 

By means of (130|) and (!3T!) the material moduli ci and C2 can be expressed in terms of 
the length £, and the constitutive moduli /z and 7 



Cl =£2(^ + ^) 



C2 



-^'(/i-7). 



(32) 



Thus, the gradation of Ci and C2 is given by the gradation of the constitutive moduli /z 
and 7. Accordingly, the governing equations for Pzx and Pzy reduce to 



[l-£2(A + 2a-V)]/3,. = /?; 
[l-f(A + 2a-V)]/3., = /?; 
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(33) 
(34) 



Due to the gradation term a -V, Eqs. ( |33|) and (IMIl may be called 'perturbed' Helmholtz 
equations. 

Differentiating equation ( 15^ with respect to y and equation with respect to x and 
using the definition of the torsion ([3]), the governing equation for the torsion (dislocation 
density) of a screw dislocation in FGM turns out to be 



[l-e{A + 2a-V)]Tzxy = T[ 





zxy 1 



where 



^ zxy 



h5{x)5{y) 



(35) 



(36) 



is the dislocation density of a straight screw dislocation in classical elasticity. Using the 
substitution 



we obtain from Eq. fl35|) the following Helmholtz equation 

h 



with 



[A - k'] ^ 



(37) 



(3^ 



(39) 



It can be seen that k is the inverse intrinsic length characteristic for the anti-plane problem 
of FGMs in the dislocation gauge approach. Therefore, 1/n is an effective length scale 
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given in terms of the gauge length scale ^ and the gradation length scale 1/a. The inverse 
length scale of gradation a is also constant. Thus, in addition to 1/a, two (constant) 
characteristic length scales, namely I and appear in the gauge-theoretical approach 
of the anti-plane problem of FGMs. The solution of (138|1 is given by 



where r = y^xM-^ and Kn is the nth order modified Bessel function of the second kind. 
Finally, the torsion of a screw dislocation in FGM using the dislocation gauge theory is 
obtained as 

r... = ^e-(-^+-^)iroM. (41) 
If we compare t he dislocation densitv (|^T1) with the corresponding one for homogeneous 



media (see, e.g.. iLazar and Anastassiadid (120091 ) ). we observe the exponential factor and 



the length scale k instead of l/£ as difference. That means the dislocation density (HTl) has 
lost the cylindrical symmetry and decays faster in the far-field in FGMs. Substituting (1411) 
into the constitutive relation f|T2l) and making use of fl30l) . we obtain for the components 
of the pseudomoment stress tensor 

Zn 

Hyz. = ~J^°^ e('^^^'+'^^^) K,{Kr) . (44) 

Due to the gradation, the fields fj4T]) - fj44|) do not possess cylindrical symmetry in contrast 
to the results for homogeneous media. It can be seen th at the expressions fHTl)-(H^ fo r 



ai — > and 02 coincide with the result obtained by iLazar and Anastassiadid (120091 ) 
for homogeneous materials. 

To solve the remaining equations (!33|) and (1M|1 . we use the constitutive relations (ITTIl 
for the components (Tzx and a^y In order to fulfill the equilibrium condition for the forces 
stresses azx,x + (^zy,y = 0, we introduce the stress function F according to 

(^zx = (/i + 'l)Pzx = -Fy , (45) 
CTzy = (yU + l)(3zy = Fx ■ (46) 

By the help of equations (j^5l) and (j46l) . we are able to express the elastic distortion (3zx 
and Pzy in terms of the stress function F according 

f3z^ = ^F^y, (47) 

Ai + 7 

Pzy = ^—F,x. (48) 

/i + 7 
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Of course, Eqs. (ST]) and (HHl) satisfy the condition (ED]). If we substitute (ST!) and (HE]) into 
(133|) and (IMl) . we find the following partial differential equation for the stress function F 

(49) 



[1 -£^(A-2a-V)]F = F°, 
where on the right hand side satisfies the equation 

[A-2a-V]F° = {iJi + -ifh5{x)5{y). 



(50) 



Eq. (150|1 is a 'perturbed' Laplace equa ti on. The solution of (l50ll . which is actually the 
Green function, reads (see iLazaii (120071 ): IWang and Pan! (120081 )) 



^0 ^ _ (/^ + 7)°fe ^(a,x+a.,)^^(^^) 

27r 

In addition, we use the ansatz for F 

and obtain from ( H9l) the Helmholtz equation for $ 

[A - K^] $ = _(^ + 7)0 65(x)(5(y) . 
The solution of (!53l) is given by 

(/i + 7)°& 



27r 



Eventually, the stress function F is calculated as 



27r 



(51) 



(52) 



(53) 



(54) 



(55) 



Substituting ( l55l) into ( l47l) and ( l48l) , we arrive at the expressions of the elastic distor- 
tions as 



In 



-{aix+a2y) ( ^ 



y — \aKi{ar) — /tJCi(/tr)] — a2[-ft'o(cf^) ~ Kq^kt 



(56) 
(57) 



The elastic distortions p.^^ and P^y are plotted in Figs. ([T]) and First of all, it can be 
seen that the elastic distortions ( 156|1 and ( l57jl are nonsingular unlike the result for FGMs 
using elasticity theory and that they are given in terms of two inverse length scales a and 
K. The component /3zj;(0,|/) has extremum values near the dislocation line at y = and 
it has a finite value at ?/ = depending on the material moduli. Unlike homogeneous 
media, it does not satisfy: PzxiS^^y) 7^ —PzxiS^i —y)- Also f3zx{x,0) possesses a maximum 
at X = and is asymmetric with respect to the plane x = 0: f3zx{x,0) ^ fizx{—^-,^\ 
Analogously, Pzy{x, 0) has extremum values near the dislocation line at x = and it has 
a finite value at a; = 0. It does not fulfill, /3zy{x,0) ^ —(3zy{—x,0), as for homogeneous 
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-10 -5 5 10 -10 -5 5 10 

(a) (b) 




(c) y/i (d) x/l 

Figure 1: Elastic distortion of a screw dislocation in FGM using gauge theory (solid line) 
and elasticity (small dashed line) are given in units of 6/[27r] and with ai = 0.1 /£ and 
a2 = 0.2/i: (a) P,^{0,y), (b) /5,^(a;,0), (c) P,y{0,y), (d) ^^,0). The dashed curves 
represent the elastic distortion in classical elasticity of a homogeneous medium. 
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Figure 3: Stress of a screw dislocation in FGM using gauge theory (solid line) and elasticity 
(small dashed line) are given in units of (yU -1-7)0 b/[2'7i] and with oi = 0.1/i and a2 = 0.2/i: 
(a) azx{0,y), (b) azx{x,0), (c) azy{0,y), (d) azy{x,0). The dashed curves represent the 
stress in classical elasticity of a homogeneous medium. 



materials. Evidently, Pzxiv, 0) has a minimum at y = and is asymmetric with respect to 
the plane y = 0: (3zy{0, y) ^ /3zy{0, —y). For FGMs the far-field of the curve of the elastic 
distortion Pzx is for negative y higher than the —1/y behavior known from homogeneous 
materials (see Fig. [T^). Due to the gradation of the material the elastic distortion in 
FGMs decays faster than in homogeneous media. 
Using (fTTj) . the stresses can be derived to be 



(Jx 



a 



_(ji^jyb ^^a,x+a,y) fV f^^^^^^) _ ^Ki{t^r)] - \Ko{ar) - Ko{tir 



yz 



271 \r 
^^~J^ ^ e("i"+''^^) (^^ [aKi{ar) - tiKi{tir)] - ai [Ko{ar) - Kq 



Kr 



(58) 
(59) 
(60) 
(61) 



The stresses azx and plotted in Figs. (jS]) and It is important to note 

that the stresses (l58l)-(l6Tl) do not possess any singularity. Due to the gradation, the 
components (!58!) - (l6Tl) have extremum values near the dislocation line and a finite value 
at the dislocation line. For FGMs the far-field of the curves of the stresses are less than 
the 1/r behavior known from homogeneous materials (see Fig. [3^ and[3li). 
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(a) (b) 

Figure 4: Three-dimensional plots of the stress of a screw dislocation in FGMs in units 
of (/X + 7)0 6/[27r] (ai = 0.1/^, aa = 0.2/£): (a) a^^ and (b) a^y. 



We note that the expressio ns (l56D-(l6Ti) coincide for ai and 02 —i- with the gauge- 
theoretical result obtained by iLazar and Anastassiadid (120091) f or homogeneous materials 
and for ai and 7 — )■ with the result given by iLazarl (120071 ) using the strain gradient 
elasticity formulation for FGMs. 



4 Conclusion 

In this paper, we have investigated a straight screw dislocation in FGMs using the gauge 
theory of dislocations. In this framework, we have derived the field equations of dislo- 
cations for general nonhomogeneous media where the constitutive moduli depend on x. 
Because of the gradation, the field equations contain gradients of the constitutive mod- 
uli. Subsequent we have specialized the field equations to exponentially graded materials. 
We found a new effective inverse length scale k characteristic for the gauge-theoretical 
anti-plane problem of FGMs. We have calculated the elastic distortion, the force stress, 
the torsion and the pseudomoment stress fields produced by a screw dislocation in FGMs 
in the framework of dislocation gauge theory. All the analytical solutions are given in 
terms of the effective length scale I/k as well as the gradation length scale 1/a. We 
found 'perturbed' Helmholtz equations as governing partial differential equations which 
we solved. We found exact analytical gauge-theoretical solutions which can be used by 
material scientists and design and manufacturing engineers. One hope is that these gauge 
solutions will provide new physical insight to possible improvements of designing FGMs. 
In the limit a — > 0, we recover well-known results of a screw dislocation in the gauge 
theory for homogeneous media. 
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